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ABSTRACT 

We present single-epoch radio afterglow observations of 24 long-duration gamma-ray burst (GRB) on a 
timescale of > 100 d after the burst. These observations trace the afterglow evolution when the blastwave has 
decelerated to mildly- or non-relativistic velocities and has roughly isotropized. We infer beaming-independent 
kinetic energies using the Sedov-Taylor self-similar solution, and find a median value for the sample of detected 
bursts of about 7 x 10^' erg, with a 90% confidence range of 1.1 x 10^" -3.3 x lO^-' erg. Both the median 
and 90% confidence range are somewhat larger than the results of multi-wavelength, multi-epoch afterglow 
modeling (including large beaming coiTections), and the distribution of beaming-corrected 7-ray energies. This 
is due to bursts in our sample with only a single-frequency observation for which we can only determine an 
upper bound on the peak of the synchrotron spectrum. This limitation leads to a wider range of allowed energies 
than for bursts with a well-measured spectral peak. Our study indicates that single-epoch centimeter-band 
observations covering the spectral peak on a timescale of St ^ 1 yr can provide a robust estimate of the total 
kinetic energy distribution with a small investment of telescope time. The substantial increase in bandwidth 
of the EVL A (up to 8 GHz simultaneously with full coverage at 1-40 GHz) will provide the opportunity to 
estimate the kinetic energy distribution of GRBs with only a few hours of data per burst. 

Subject headings: gamma-rays:bursts 



1. introduction 

The energy budget of gamma-ray bursts (GRBs) provides 
fundamental insight into the nature of the explosions, the re- 
sulting ejecta properties, and the identity of the central com- 
pact remnant ("engine"). While the isotropic -equivalent 7- 
ray energy (fi-y.iso) can be easily determined from a mea- 
surement of the burst fluence and redshift, a complete ac- 
counting of the energy budget requires detailed observations 
of the afterglow emission. The afterglow observations pro- 
vide a measure of the isotropic-equivalent blastwave kinetic 
energy (iiK.iso), as well as the explosion geometry (quanti- 
fied by a jet opening angle, 6j). The resulting beaming cor- 
rections, = 1 - cos{dj), can be substan tial, approachin g 
three orders of magnitude in some cases (Frail etal."2001f 
Panaitescu & Ku mar 2002; Be rger et al . 2003a; Bloom et aT| 
2003). To properly determine ^K.iso and fi, it is essential to 
observe the afterglows from radio to X-rays over timescales 
of hours to weeks, clearly a challenging task. This is particu- 
larly a problem for the subset of "dark" GRBs for which the 
lack of detected optical emission, or large extinction, prevent 
a determination of /iK.iso and likely fi, (e.g.. lBerger et al.l2002l; 
IPiro et al. 2002). 

Over the past decade detailed afterglow observations 
have been obtained at a great cost of telescope time for 
about 20 long-duration GRBs, with the basic result that 
the beaming coiTections are large and diverse, leading to 
typical true energies of E-., ^ Ek '--^ 10^' erg j Frajl et alj 
2001l:lPanait escu & Kumar 2001, 2 002HBerger"eta l. 2003a b; 
Bloom et al.l ,200 3). More recently, it has been recognized 



that some nearby long GRBs have much lower energies, 
Ejso ~ 10"^^ - 1 0^" erg, and appear to be quasi-isotropic 
(iKulkarni et a l.l[l998; Soderberg et al. 2004b, 2006). Simi- 
larly, some bursts appear to have large beaming-corrected en- 

' Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cam- 
bridge, MA 02138 



ergies of ~ 10^^ erg (ICenko et alJl2010"bllah . The existence 
of these highly energetic bursts depends at least in part on 
the ability to correctly infer their large beaming corrections. 
Indeed, the inference of jet opening angles from breaks in 
the afterglow light curves has become controversial in re- 
cent ye ars due to conflicting trends in optical a nd X-ray light 
curves (iLiang et aLllIOOSl: iRacusin et alJl2009l) . Similarly, in 
some cases a two-component jet has been inferred, with a nar- 
row core dominating the 7-ray emission and a wider compo- 
nent dominating the afterglow emission ( iBerger et alj|2003bt 
iRacusin et al.ll2008h . Numerical simulations suggest that off- 
axis viewing angles can also le ad to shallow b reaks that may 
be missed or mis-interpreted (van E erten et al.ll2010 ). 

In addition to potential difficulties with the inference of 
the 7-ray and kinetic energies measured from the early after- 
glow emission only pertain to the relativistic ejecta. The exis- 
tence of a substantial component of mildly relativistic ejecta 
can only be determined from observations at late times when 
such putative material can refresh the forward shock. Clearly, 
the existence of substantial energy in a slow ejecta compo- 
nent will place crucial constraints on the activity lifetime of 
the central engine. 

Such late-time observations also have the added advantage 
that they probe the blastwave when it has decelerated to non- 
relativistic velocities and hence rou ghly approaches isotropy 
dFrail et al.ll2000t iLivio & Waxmanli200 0). This aflows us to 
use the well-established Sedov-Taylor self-similar solution, 
with negligible beaming corrections, to estimate the total ki- 
netic energy of both the decelerated ejecta and any additional 
initially non-relativistic material. Since the peak of the af- 
terglow spectrum on these timescales is located in the radio 
band, the lack of optical afterglow emission (e.g., due to ex- 
tinction) does not have an effect on the ability to determine 
Ek- , , , , 

This approach was first exploited by Frail et al. (2000) to 
model the late-time radio afterglow emission of GRB 970508 
(at (5f > 100 d) from which the kinetic energy was infeiTed to 
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he Ek ^ 5x 10^° erg. iBerger et all ( |2004|) used the same ap- 
proach to model the radio afterglow emission of GRB 980703 
on timescales of > 40 d, and to re-model GRB 970508. They 
found k inetic energies of Eg ^ 3 x 10^^ erg for both bursts. 
Finally, iFrail et alT (12005) modeled the radio emission from 
GRB 030329 at <5f > 50 d and found Ek - 10^' erg. Only 3 
bursts have been studied in this fashion so far because only 
those events have well-sampled radio light curves on the rele- 
vant timescales of St > 100 d. 

However, the kinetic energy can still be estimated using the 
same methodology even from fragmentary late-time radio ob- 
servations. Such an approach will naturally result in larger 
uncertainties for each burst, but it can be applied to a much 
larger sample of events. Here we present such an analysis for 
24 long-duration GRBs with radio observations at > 100 d, 
but with only 1-3 data points (at 1.4 to 8.5 GHz) per burst. 
Using these observations we infer robust ranges for the ki- 
netic energy of each burst and for the population as a whole. 
The plan of the paper is as follows. The radio observations 
are summarized in ^ The model for synchrotron emission 
from a Sedov-Taylor blastwave, and the various assumptions 
we employ are presented in ^ In ^we detail the resulting 
kinetic energies and the range for the overall sample, and we 
compare these results to multi-wavelength analyses of early 
afterglows in ^ We conclude with a discussion of future 
prospects. 

2. RADIO DATA 

We use radio observations of 24 long GRBs at St > 100 
d since on those timescales the blastwave is expected to be- 
come non-relativistic and roughly isotropic (i Livio & WaxmanI 
120001) . and the peak of the afterglow emission is at or be- 
low the centimeter band. This has been confirmed with de- 
tailed data in the case of GRBs 970508, 980703 and 30329 
(IFrail et al.1 120001: IBerger et al.1 l200l [prail et alJ l2005h . We 
restrict the analysis to GRBs with a known redshift and with 
early-time detections, which for the case of a single detection 
or upper limit allow us to infer that the peak of the spectrum 
has transitioned below our observing frequency. 

The observations are primarily from the Very Large Ar- 
ray (VLA^), with the exception of GRBs 980425 and 011121 
which were observed with the Australia Telescope Com- 
pact Array (ATCA). The data were obtained between 1997 
and 2009 as part of a long-term GRB radio program (e.g., 
IFrail et al.ll2003ah . 

For the purpose of our analysis, we separate the bursts into 
three categories based on the quality of the data. In Group A 
are 3 bursts with late-time detections at multiple frequencies 
that constrain the peak of the synchrotron spectrum (the same 
three bursts that have been studied in d etail by IFrail et alj 
l2000HBerger etani200llF^ail et al.ll2005h . In Group B are 1 1 
bursts with single-frequency detections, while Group C con- 
sists of 10 GRBs with late-time non-detections. The VLA 
measurements and relevant burst properties are listed in ta- 
ble[I] 

3. SYNCHROTRON EMISSION FROM A NON-RELATIVISTIC 
BLASTWAVE 

Our m odeling of t he ra dio d ata follows t h e met hod- 
ology of IFrail et alj (l2000h and IBerger et alj (|2004|) for 

^ The National Radio Astronomy Observatory is a facility of the National 
Science Foundation operated under cooperative agreement by Associated 
Universities, Inc. 



the case of a uniform density medium . For the typi- 
cal expected parameters of long GRBs, the initially coUi- 
mated blastwave approaches spherical symmetry and decel- 
erates to non-relativistic velocity on similar timescales, w 

150(£K.iso.52 /«.)^/'f- d and t ^R « 40(£K,iso,52/«.)'/'f- d, 
respectively dLivio & Waxmanll2000l) : here, n,, is the circum- 
burst density in units of cm"^ and tj is the "jet break" time 
at which the jet begins to expand sideways (i.e., r(fj) ^ 9^^, 
where F is the bulk Lorentz factor). In this paper we as- 
sume that the blastwave has transitioned to the non-relativistic 
isotropic phase by the time of our observations and subse- 
quently check for self-consistency. 

The blastwave dynamics in the non-relativistic phase are 
described by the Sedov-Taylor self-similar solution with 
r(t) oc {Esjt^/ny/^. To calculate the synchrotron emission 
emerging from the shock-heated material, we make the usual 
assumptions: (i) the electrons are accelerated to a power- 
law energy distribution, A^(7) oc 7"'' for 7 > 7,,,, where 7^ 
is the minimum Lorentz factor; (ii) the value of p is 2.2 as 
inferred from several b ursts (e.g.. lPanaitescu & KumaHr2001l 
I2OO2I: lYost et al.1 [20031) : and (iii) the energy densities in the 
magnetic field and electrons are constant fractions (eg and 
respectively) of the shock energy density. Accounting for 
synchrotron emissivity and self-absorption, and including the 
appropriate redshift transformatio ns, the flux observed at fre- 
quenc y v and time t is given bv (IFrail et al.ll2000l: IBerger et alj 
12001 : 

F, = Fo(t/tQrn(l+z)iyfHl-e-^-^)Miy/iy„df2\'y/'ym), (D 
where the optical depth is given by: 

= To(f /fo)"- [(1 +z>]-'-P^^^^f2(iy/iy,n), (2) 

the synchrotron peak frequency, corresponding to electrons 
with 7 = 7,„, is given by: 

i^,„ = iyo(t/tor'-(l+zT\ (3) 
and the function //(x) is given by 

Mx)= f F{y)/P-'^^^dy, (4) 
Jo 

where F(y) is an inte gration over Bessel functions 
dRvbicki & Lightmanll9 79). The temporal indices in the case 
of a uniform density medium are af = 11/10, Ofr = 1-3/9/2, 
and am = -3. The normalizations are such that Fq and tq are 
the flux density and optical depth at a frequency of 1/ = 1 Hz 
at f = ?(), and I'o is the synchrotron peak frequency in the rest 
frame of burst at f = to. Furthermore, the synchrotron self- 
absorption frequency, Va, is defined by the condition T^(i'fl) = 
1. 

We fit this synchrotron model to our radio data using Fq, 
Tq, and 1^0 as free parameters. Since we have no a priori 
knowledge about the expected values of the synchrotron spec- 
trum parameters we assume that they follow a flat distribution 
in log-space. We note that any further assumption about the 
distribution of these parameters will only serve to restrict the 
resulting energy distributions, and we therefore consider our 
assumed flat distribution to be conservative. For the detected 
objects (Groups A and B) we retain all solutions that repro- 
duce the measured flux density within the error bars, while 

' Since we use a single epoch of observations for each GRB our inferred 
density can be easily converted to a mass loss rate for the case of a wind 
medium. The difference in dynamical evolution between these two models 
does not have an effect in this case. 
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for the non-detections we use 3(t as an upper bound. More- 
over, for the bursts in Group A, z/q is constrained by the muhi- 
frequency observations and no additional constraints are re- 
quired. However, for the bursts in Groups B and C, which 
have only a single-frequency observation, we require that both 
Vm and Va have values below the observing frequency since the 
light curves are always declining at the time of our observa- 
tions'*. 

Using the allowed ranges of Fq, tq, and V(, we determine 
the set of relevant physical parameters: n^, 7,,,, and B, where 
B is the magnetic field strength . The radius of the blastwave, 
r, rem ains unconstrained (e.g., iFrail et al]l200(]l: iBerger et"aD 
I2OOI : 

B = nj{p + 2r%^_,2(ru/dL2i)' G, (5) 

7„, = 6.7(/:> + 2)Fo-52i'o,i^('-i7/'iL,28r^, (6) 

He = 3.6 X l0^%,r]iFQ_^2'^Qg''''^^T032ri](rn /dL,2&)~^ cm"^ 

(7) 

c„ = (1.67 X 10V(5.4 X 102)(i-''>/2(/, + 2)V(/7- 1) (8) 

where £^^ 28 is the luminosity distance in units of 10^^ cm as- 
suming the standard cosmological parameters (Ho = 71 km s"' 
Mpc"\ = 0.27, and ft^ = 0.73), and ry is the reciprocal of 
the thickness of the emitting shell. 

The unknown radius of the blastwave can be con- 
strained by introducing a relationship between Est = 
ngmp(r/l.Q5)^[ti^i{/(l+z)]~~ and the energy in the electrons 
and the magnetic field. We use the condition that at most 
half of the blastwave energy is available for accelerating 
electrons and producing the magnetic field, i.e., (EB + Et,) < 
Est/ 2. The total energy in the accelerated electrons is E,, = 
[(/?- 1 )/(/?- 2)]«j,7,„c^y, while the energy in the magnetic 
field is Eb = -S^V/Stt, where V = Aur^ /rj is the volume of 
the synchrotron emitting shell. The energy budget is min- 
imized near equipartition (i.e., E^ ~ Eb), and we use this 
constraint to de t ermine the minimu m required energy (e.g., 
iFrail et al.l l2000'.'Berger et al."2004'): this conclusion was ver- 
ified with radio interferometric measurements of the size of 
GRB 030329 ( Tavlor et al.il2004t IFrail et al.li2005i) . 

Finally, using the inferred radius for each possible solu- 
tion, we require for self-consistency that /3 J; 1, where /3 = 
2r(l +z)/5ct. The resulting (3 distributions are shown in Fig- 
ure 12] These results indicate that most of the detected bursts 
obey the self-consistency requirement, although we reject 
GRBs 000926, 020819, and 021004 from the analysis since 
> 50% of their allowed solutions lead to relativistic velocities 
at the time of the observations. This does not rule out that 
the Sedov-Taylor solution is applicable, but simply indicates 
that additional observations are required to narrow down the 
range of allowed solutions. We furthermore find that the bulk 
of the upper limits do not rule out relativistic expansion, and 
we therefore do not use these limits in the energy distribution 
analysis below. 

4. THE DISTRIBUTION OF GRB KINETIC ENERGIES 

The resulting solutions for each burst can be cast in terms 
of a two-dimensional parameter space in Ug versus Est- Thus, 
there is a degeneracy between the two parameters, in the sense 
that larger densities lead to lower energies. Clearly, the bursts 
in Group A, for which the peak of the synchrotron spectrum 

* The opposite case of either or Va being larger than the observing 
frequency leads to rising light curves. 



is well-determined, lead to the best constraints in this two- 
dimensional phase space. Indeed, as shown for GRB 980703 
(Figure [Til, the allowed range of energies for all solutions 
that reproduce the observed flux density is about lO'*^- lO^-' 
erg, while the solutions that also satisfy the requirements 
that (EB + Ee) < Esj/2 span a much narrower range of about 
lO^i - 10^2 gj.g^ ^jth 

a roughly log-normal distribution cen- 
tered on log(£'sT) ~ 5 1 .6. 

The bursts with only single-frequency observations (Group 
B) cover a much larger area in the Esj - phase-space since 
only an upper bound can be placed on v„, and Va- To fur- 
ther constrain the energy we place an additional conserva- 
tive constraint on the density of rie < 100 cm"^, motivated 
by the results of detailed broad-band modeling tha t shovy 
n, ~ 0.1- 10 c m"^ (e.g., Panaite scu & Kumad I200TI 120021: 
I Yost et ani2003h . Given the anti-correlation between density 
and energy, our conservative limits lead to a wider range of 
allowed energies than if we chose a limit of «(. < 10 cm""*. 
An example of this additional constraint for a Group B burst 
(GRB 010921) is shown in Figure [T] We do not place a lower 
bound on the density, since for the phase-space of allowed 
solutions this would not lead to a significant change in the 
energy distribution. We stress that beyond placing an upper 
bound on no constraints have been placed on the distribu- 
tions of either n,, or (3 since both are inferred, and not input, 
parameters in our model. 

As noted above, we do not consider the energies for the 
bursts in Group C since the upper limits generally allow a 
wide range of solutions that are not consistent with the Sedov- 
Taylor formulation. This indicates that the limits are generally 
not deep enough to provide a meaningful constraint on the 
energy. Future deep radio observations may provide much 
better constraints (see below). 

The resulting energy probability distributions for the bursts 
in Group A and Group B are shown in Figure [3] The median 
energy and 90% confidence range (i.e., 5-95% of the dis- 
tribution) for each burst are listed in Table |2] We include in 
these ranges the small subset of solutions that lead to (3 values 
in slight excess of 1 since these are at most mildly relativistic 
and furthermore do not significantly change the distributions 
(Figure |3]l. We find that varying the electron power law index 
over the range /? = 2.1-2.5 (e.g. JCurran et al.l2008l) leads to a 
change in the median energy of only 0.1-0.2 dex (compared 
to our fiducial value of p = 2.2), with larger values of p lead- 
ing to lower median energies. Similarly, varying the magnetic 
energy fraction away from equipartition to = 0.1 and 0.01, 
leads to an increae in the median energy of about 0.25 and 
0.5 dex, respectively. Both of these effects are much smaller 
than the overall spread in energy for each burst, but they do 
produce minor systematic trends. 

The combined distribution for the subset of 1 1 bursts whose 
solutions are generally self-consistent is shown in Figure |4] 
The median and 90% confidence ranges are 7 x 10^' erg and 
1.1 X 105'*-3.3 X lO^^^erg. 

5. DISCUSSION AND CONCLUSIONS 

The key results of our analysis are that the median en- 
ergy for the 1 1 bursts with self-consistent solutions is Ek ~ 
7 X 10^' erg, while the 90% confidence range is 1.1 x 10^°- 
3.3 X lO^^* erg. The median value is about a factor of 3 
times higher than previous calorimetric measurements for 
GRBs 970508, 980703, and 030329, for which energies of 
3 X 10^', 3 X 10^', and 10^' erg, respe ctively, were determined 
(IBerger et al.ll2004t IFrail et al.ll2005l) . 
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Similarly, the inferred energies are somewhat larger than 
the distributions of beaming-corrected 7-ray and kinetic en- 
ergies inferred from broad-band multi-epoch studies (Fig- 
ure |5]l. From various su ch analyses, the median 7-ray energ 
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is » 8 X 10^° erg (iFrail et al.l '2001; Bl oom et alJ 12003 
iFriedman & Bloomll2005h . while the median kinetic energy is 
(Ek) ~5 X 10^0 erg (e.g.. iPanaitescu & Kumadl200lL 120021: 
lYost et alJl2003l) : see Figure |5] In both cases the 90% range 
spans about 2.5 orders of magnitude, somewhat narrower than 
our inferred 90% confidence range for Est- The extension 
to larger energies found in our analysis mainly reflects the 
lack of spectral peak determinations for the bursts with single- 
frequency observations (see Figure O. These large energies 
can be generally eliminated with a measurement of the syn- 
chrotron peak in the GHz frequency range (e.g.. Group A 
bursts; Figure |5]). 

In the con text of our results we note that recent numeri- 
cal work by IZhang & MacFadvenI ( 1200 9) led these authors 
to conclude that the timescale to reach isotropy is ^ 10^ yr 
rat her than ^1 yr a s indicated by the analytic formulation 
of iLivio & Waxmanl (2000) which we follow here. As a re- 
sult, they note that using the Sedov-Taylor formulation may 
lead to an erroneous estimate of the kinetic energy. However, 
inspection of the resulting potential disrepancies reveals that 
this effect is at most a factor of 2 as long as self-consistency 
between the inferred energy and density and the transition 
to the Sedov-Taylor phase is ensured (see their Figure 10). 
The discrepancies become larger if the wrong timescale is as- 
sumed for the transition to non-relativistic expansion, but this 
quantity is not a free parameter. Indeed, our distributions of /? 
values point to self-consistency for most bursts, and allow us 
to reject objects that are potentially still relativistic. Since the 
potential systematic uncertainty of about a factor of 2 is sig- 
nificantly smaller than the overall spread in allowed energy 
for each burst, we do not consider this to be an obstacle to our 
analysis, or to future work on the energy scale using late-time 
radio measurements. 

As clearly demonstrated in Figure [3] the most constrained 
energy determinations require a measurement of the syn- 
chrotron spectral peak (Group A); the absence of such a con- 



straint requires additional assumptions about the circumburst 
density and results in a much wider energy range. Indeed, 
this is the key reason for the wider range of allowed high en- 
ergy solutions (> 10^^ erg) compared to the results for E^ 
and Ek (Figure |5]l. Observations of GRBs 970508, 980703, 
and 030329 demonstrate that the spectral peak is typically lo- 
cated at ^ few GHz on a timescale of 150 d. Thus, ob- 
servations in the 1-10 GHz range on a timescale of ^ few 
hundred days should allow us to determine the peak flux and 
frequency. This will in turn provide an energy estimate with 
a similar level of precision to the results of early-time broad- 
band modeling. 

This is a fortuitous conclusion since with the full frequency 
coverage of the Expanded VLA (EVLA) it will soon be pos- 
sible to cover this entire range in a few hours of observations 
to a sensitivity that is about an order of magnitude better than 
the VLA. As we demonstrated here, such a modest invest- 
ment of observing time (2-3 hours per burst) can yield a ro- 
bust estimate of the GRB energy distribution, regardless of 
the ability to measure jet opening angles. Pursuing these ob- 
servations for all bursts with a measured redshift will require 
only ^50-100 hr of EVLA time per year Indeed, with such 
observations we should be able to constrain the energy dis- 
tribution to a comparable level as existing studies within a 
single year given that about 30 GRBs with known redshifts 
just from 2009 are now available for EVLA observations (a 
similar sample is available from 2008 bursts). In the longer 
term, the large number of objects will allow us to test the en- 
ergy distribution as a function of redshift, at least over the 
range z ^ 1 - 3 where the bulk of the detected bursts occur 
(Berge ret al.l l2005; Jakobsson et al. 2006). Similarly, this ap- 
proach will be particularly useful for bursts that lack detailed 
optical or X-ray light curves due to observational constraints 
or dust extinction, and for bursts with controversial estimates 
of the jet opening angles. 



We thank Dale Frail and Eli Waxman for helpful discus- 
sions and comments on the manuscript. 
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TABLE 1 

Late-time Radio Afterglow Measurements 



GRB 


z 


St 
(d) 


V 

(GHz) 


(MJy) 


Ref. 




0.835 


1 17.55 


8.46 


355 ± 47 


rrail et al. (2000) 






1 17.55 


4.86 


425 ± 57 






117.55 


1.43 


206 ± 63 




970828 


0.958 


157.99 


8.46 


< 51 


Diorsovski et al. (2001) 


980425 


0.0085 


248.20 


8.70 


700 ± 200 


Kulkami et al. (1998") 


980703 


0.966 


143.79 


8.46 


110 ±20 


hrail et al. (2003b) 






143.79 


4.86 


146 ±24 








134.85 


1.43 


99 ±25 




990506 


1.307 


141.23 


8.46 


< 75 


Taylor et al. (2000) 


991208 


0.706 


291.58 


8.46 


51 ± 15 


Galama et al. (2003^ 


000210 


0.846 


108.37 


8.46 


< 78 


_Frail et al. (2003a) 


00030 IC 


2.030 


506.10 


8.46 


39± 11 


^rger et al. (2000) 


000418 


1.118 


405.76 


8.46 


38± 11 


^erger et al. (2001) 


ouoyi 1 


1 ACQ 


125.78 


0.40 


< 54 


Frice et al. (ZUUz) 


000926 


2.066 


257.43 


8.46 


75 ± 21 


garrison et al. (zUOlj 


010222 


1.477 


206.63 


8.46 


< 42 


Frail et al. (2003a) 


010921 


0.451 


225.42 


8.46 


52 ± 15 


rrail et al. (2003a) 


011121 


0.362 


132.07 


8.70 


< 141 


Frail et al. (2003a) 


none 1 Q 
uzuoiy 


n /1 1 1 




O.40 


ly ± Zj 


ToVi-ii-^ccz-in at ni ^^^^^^^^ 
jaKODSson et ai. [ zuu j ) 


021004 


2.329 


140.24 


8.46 


94 ±16 


Frail et al. (2003a) 


030226 


1.986 


113.85 


1.43 


< 117 


Frail et al. (20g3a) 


030329 


0.168 


135.48 


8.46 


1525 ±56 


Frail et al. (2005) 






129.57 


4.86 


1955 ±62 






129.58 


1.43 


1276 ±56 




031203 


0.105 


137.15 


8.46 


426 ±37 


Soderberg et al. (2004b) 


0504 16A 


0.654 


182.28 


8.46 


< 114 


Soderberg et al. (2006) 


070125 


1.547 


341.96 


8.46 


64±18 


Chandra et al. (2008) 


0706 12A 


0.617 


488.54 


8.46 


101±39 


Frail et al. (2003a) 


090323 


3.570 


131.18 


8.46 


<81 


Cenko et al. (2010a) 


090902B 


1.822 


199.16 


8.46 


<48 


Cenko et al. (2010a) 



Note. — " Limits are 3cr. 
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TABLE 2 

GRB Energies Inferred from Calorimetry 



GRB 


(log(£sT)> 


log(£sT).90'' 


log(£^)* 






(erg) 


(erg) 


(erg) 


(erg) 


Group A 


970508 


51.8 


51.3-52.5 


50.6 


51.3 


980703 


51.6 


51.1-52.2 


51.0 


51.5 


030329 


51.3 


50.9-51.8 


49.9 


50.4 


Group B 


980425 


49.4 


48.9-50.1 


47.8 


~50 


991208'' 


51.9 


51.0-53.4 


51.2 


50.4 




52.1 


51.1-53.6 






000301C 


52.4 


51.8-53.5 


50.9 


50.5 




52.6 


51.8-53.9 






000418 


50.6 


49.8-51.6 


51.7 


51.5 


000926 


52.0 


51.6-52.8 


51.2 


51.2 




52.7 


51.8-54.3 






010921 


51.6 


50.6-53.1 


<51.2 






51.9 


50.7-53.7 






020819 


51.3 


50.4-52.6 


< 51.8 






51.8 


50.6-53.6 






021004 


51.8 


51.4-52.3 


50.9 






52.8 


51.7-54.5 






031203 


51.1 


50.2-52.5 


49.5 


49.2 




51.5 


50.3-52.5 






070125 


52.2 


51.6-53.2 


52.4 


51.2 




52.6 


51.7-54.0 






0706 12A 


52.1 


51.6-53.2 


<52.0 





Note. — " This is the 90% confidence range for the energy of each 

burst. 

^ Values for E-^ are taken f rom lFrail et alj |2001|) . |Bloom et al] 120031) . and 
iFried man & Bloom' i'200j. 

^ Val ues for Ek are tak en fr om[Panaitescu_&Kumaj (200^, [Bergeretal 
2003b) . lYostetaL aOOl. iSoderberg et alJ I2004H) . iSoderberg et al. 
2w3i . ISoderberg etalJ I2006^ . iCenko et al] 1201 OH) , and iCenko etal 

^ The first line is for a strict cut-off of /3 < 1, while the second line allows 
a small fraction of solution with /3 slightly larger than 1 (Figure [3]. 
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48 



49 



50 



51 



52 



53 



logfEgj] (erg) 



48 



49 



50 



51 



52 



53 



logfEgj] (erg) 



Fig. 1 . — Electron number density plotted against kinetic energy for two representative cases. The light gray regions indicate the phase-space that leads 
to a predicted flux density in agreement with the observed values. The medium gray regions encompass the subset of solutions that satisfy the condition 
(Es+Ee) < Est/2. The black regions marks the subset of solutions that satisfy /3 < 1 in the Sedov-Taylor framework. Left: Group A burst with a well-defined 
spectral peak. Right: Group B burst with a single frequency detection for which we use the additional limit that n < 100 cm~^. This figure highlights the 
significant advantage of measuring the spectral peak. 
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Fig. 2. — Normalized histograms of inferred expansion velocity, f) = v/c, at the time of our observations. A value of < 1 (vertical hnes) is required for 
self-consistency, and this is indeed the case for the bulk of the acceptable solutions. Note that the scales for the three groups are different. 
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Fig. 3. — Normalized histograms of GRB energies calculated using the Sedov-Taylor solution. The light gray histograms for some Group B bursts indicate the 
subset of solutions with a strict cut-off of /3 < 1. GRBs 000926, 020819, and 021004 are rejected from our sample since the bulk of their solutions still lead to 
relativistic expansion. 
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0.15 




log[E] (erg) 

Fig. 4. — Normalized distribution of GRB kinetic energies calculated using the Sedov-Taylor solution for the sample of 11 bursts with self-consistent solutions 
(Figure[3). The median and 90% confidence range are marked by a horizontal bar. 
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log[E] (erg) 

Fig. 5. — Normalized distributions of GRB kinetic energies calculated using the Sedov-Taylor solution (black), and for the subset of 3 bursts in Group A 
(gray ). Also shown for comparison are the distributions of beaming-coirected 7-ray energies (red: hatch = known 6j values; o pen = 6, lower or upper limits; 
IRriedman & Bloom 2005) and beaming-corrected kinetic energies from broad-band early afterglow modeling (blue; Panaite scu & Kumaj|2002l ; IBerger et al.l 
|2003b; Yost et al. 2003; Soderberg et al. 2004b a, 2006; Cenko et al. 2010b a). The median and 90% confidence range for each energy component are marked 
by a horizontal bar Our inferred median energy and 90% confidence range are larger than the median of both Ej and Ek, but this is mainly due to the bursts in 
Group B for which the spectral peak is not measured. Future observations with the EVLA will lead to much tighter constraints (see gray histogram) for a larger 
sample. 



